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As neutrino long baseline experiments enter a new domain of precision, the careful study 
of systematic errors due to poor knowledge of production cross sections for pions and 
kaons require more dedicated measurements for precise neutrino flux predictions. The 
cosmic ray experiments require dedicated hadron production measurements to tune sim- 
ulation models used to describe air shower profiles. Among other goals, the NA61/SHINE 
(SPS Heavy Ion and Neutrino Experiment) experiment at the CERN SPS aims at preci- 
sion measurements (5% and below) for both neutrino and cosmic ray experiments: those 
will improve the prediction of the neutrino flux for the T2K experiment at J-PARC and 
the prediction of muon production in the propagation of air showers for the Auger and 
KASCADE experiments. NA61/SHINE took data during a pilot run in 2007 and in 2009 
and 2010 with different carbon targets. The NA61/SHINE set-up and spectra for positive 
and negative pions obtained with the 2007 thin (4% interaction length) carbon target 
data are presente dP 

Keywords: p+C interaction, inelastic cross section, inclusive pion spectra 



1. Physics motivation 

The NA61/SHINE (SPS Heavy Ion and Neutrino Experiment) experiment at the 
CERN SPS pursues a rich physics program in various fields .I^HslFii-st^ precise hadron 
production measurements are performed for improving calculations of the neutrino 
flux in the T2K neutrino oscillation experiment,'^ as well as for more reliable simula- 
tions of cosmic-ray air showers in the Pierre Auger and KASCADE experimentsEEl 
Second, p+p, p+Pb and nucleus-l-nucleus collisions will be studied extensively at 
SPS energies. This article presents first NA61/SHINE results on charged pion spec- 
tra in p-|-C interactions at 31 GeV/c which are needed for an accurate neutrino 
flux prediction in the T2K experiment. The results are based on the data collected 
during the first NA61/SHINE run in 2007. 

T2K is a long baseline neutrino experiment in Japan, which uses a high inten- 
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sity neutrino beam produced at J-PARCQ It aims to precisely measure the — )■ i/g 
appearance and disappearanceP In order to generate the neutrino beam a high 
intensity 30 GeV (kinetic energy) proton beam impinging on a 90 cm long graphite 
target is used, where tt and K mesons decaying into (anti)neutrinos are produced. 
The neutrino fluxes and spectra are then measured both at the near detector com- 
plex, 280 m from the target, and by the Super-Kamiokande (SK) detector located 
295 km away from the neutrino source and 2.5 degrees off-axis. Neutrino oscillations 
are probed by comparing the neutrino flux measured at SK to the predicted one. 
In order to predict the flux at SK one uses the near detector measurements and ex- 
trapolates them to SK with the help of Monte Carlo simulations. Up to now, these 
Monte Carlo predictions are based on hadron production models only. For more pre- 
cise predictions, which would allow the reduction of systematic uncertainties to the 
level needed for the T2K physics goals, measurements of pion and kaon production 
off carbon targets are essential.i^: The purpose of the NA61/SHINE measurements 
for T2K is to provide this information at exactly the proton extraction energy of 
the J-PARC Main Ring synchrotron, namely 30 GeV kinetic energy (approximately 
31 GeV/c momentum). Presently, the T2K neutrino beam-line is set up to focus 
positively charged hadrons, in such a way that it produces a i^^ beam. Spectra 
of positively charged pions presented in this paper constitute directly an essential 
ingredient in the neutrino flux calculation. 

2. The NA61/SHINE detector 

The NA61/SHINE experiment is a large acceptance hadron spectrometer in the 
North Area H2 beam-line of the CERN SPS. The schematic layout is shown in 
Fig. [l] together with the overall dimensions. 

The main components of the current detector were constructed and used by the 
NA49 collaboration.'^ A set of scintillation and Cherenkov counters as well as beam 
position detectors (BPDs) upstream of the spectrometer provide timing reference, 
identification and position measurements of the incoming beam particles. The main 
tracking devices of the spectrometer are large volume Time Projection Chambers 
(TPCs). Two of them, the vertex TPCs (VTPC-1 and VTPC-2 in Fig. [l]), are 
located in a free gap of 100 cm between the upper and lower coils of the two super- 
conducting dipole magnets. Their maximum combined bending power is 9 Tm. In 
order to optimize the acceptance of the detector at 31 GeV/c beam momentum, the 
magnetic field used during the 2007 data taking period was set to a bending power of 
1.14 Tm. Two large TPCs (MTPC-L and MTPC-R) are positioned downstream of 
the magnets symmetrically to the beam line. The TPCs are filled with Ar:C02 gas 
mixtures in proportions 90:10 for VTPCs and 95:5 for MTPCs. The particle identi- 
fication capability of the TPCs based on measurements of the specific energy loss, 
dE/dx, is augmented by time-of-flight measurements using Time-of-Flight (ToF) 
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Fig. 1. (Color online) The layout of the NA61/SHINE experiment at the CERN SPS (top view, 
not to scale). The chosen right-handed coordinate system is shown on the plot. The incoming 
beam direction is along the z axis. The magnetic field bends charged particle trajectories in the 
X — z (horizontal) plane. The drift direction in the TPCs is along the y (vertical) axis. 

detectors. The ToF-L and ToF-R arrays of scintillator pixels have a time resolu- 
tion of better than 90 ps.f^ Before the 2007 run the experiment was upgraded with 
a new forward time-of-flight detector (ToF-F) in order to extend the acceptance. 
The ToF-F consists of 64 scintillator bars with photomultiplier (PMT) readout at 
both ends resulting in a time resolution of about 115 ps. The target under study is 
installed 80 cm in front of the VTPC-1. The results presented here were obtained 
with an isotropic graphite target of dimensions 2.5(W)x 2.5(H) x 2 (L) cm and with 
a density of p = 1.84 g/cm'^. The target thickness along the beam is equivalent to 
about 4% of a nuclear interaction length (Aj). 

3. Analysis techniques 

This section presents the procedures used for data analysis to extract pion cross 
sections. Crucial for this analysis is the identification of the produced pious. De- 
pending on the momentum interval, different approaches have been adopted, which 
lead also to different track selection criteria. The calibrated dE/dx distributions as 
a function of particle momentum for positively and negatively charged particles are 
presented in Fig. [2] 

The task is facilitated for the negatively charged pions, by the observation that 
more than 90% of primary negatively charged particles produced in p+C interac- 
tions at this energy are 7r~, and thus the analysis of tt^ spectra can also be carried 
out without additional particle identification. 

In the low momentum region (less than about 1 GeV/c), it is sufficient to dis- 
tinguish pions from electrons/positrons, kaons and protons by means of particle 
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Fig. 2. (Color online) Specific energy loss in the TPCs for positively (right) and negatively (left) 
charged particles as a function of momentum. Curves show parameterizations of the mean dE/dx 
calculated for different particle species. 



identification via measurements of specific energy loss {dE/dx) in tfie TPCs. A 
reliable identification of 7r+ mesons was not possible at momenta above 1 GeV/c 
where the Bethe-Bloch (BB) curves for pions, kaons and protons cross each other. 
On the other hand, for n~ mesons, where the contribution of K~ and antiprotons 
is almost negligible, the dE/dx analysis could be extended in momentum up to 
3 GeV/c allowing consistency checks with the other analysis methods in the region 
of overlap. 




p [GeV/c] dE/dx [mip] 



Fig. 3. (Color online) Right: Mass squared, derived from the ToF-F measurement and the fitted 
path length and momentum, versus momentum p. The lines show the expected mass squared values 
for different particles. Left: Example of two-dimensional m?-dE / dx plots for positively charged 
particles in the momentum range 2-3 GeV/c. 2a contour around fitted pion peaks are shown. 



High purity particle identification can be performed by combining the to/ and 
dE/dx information. Moreover, in the momentum range 1-4 GeV/c, where dE/dx 
bands for different particle species overlap, particle identification is in general only 
possible using the tof method, see Figs. [3j The ToF-F detector was designed to 
cover the necessary acceptance in momentum and polar angle required by the T2K 
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experiment, although hmited to particle momenta above about 0.8 GeV/c. 
Indeed, three analysis methods were applied to obtain pion spectra: 

(1) analysis of 7r~ mesons via measurements of negatively charged particles {h~ 
analysisf^ 

(2) analysis of 7r+ and tt^ mesons identified via dE/dx measurements in the TPCs 
{dE/dx analysis at low momentum^ and 

(3) analysis of tt"*" and ■n~ mesons identified via time-of- flight and dE/dx measure- 
ments in the ToF-F and TPCs, respectively (to/ — dE/dx analysis)^ 

Each analysis yields fully corrected pion spectra with independently calculated 
statistical and systematic errors. The spectra were compared in overlapping phase- 
space domains to check their consistency. Complementary domains were combined 
to reach maximum acceptance. 

The agreement between the spectra obtained by different methods is, in general, 
better than 10%. Note, that data points in the same {p, 6) bin from different analysis 
methods are statistically correlated as they result from the analysis of the same data 
set. In order to obtain the final spectra consisting of statistically uncorrelated points 
the measurement with the smallest total error was selected. 

The analysis was based on a sample of 521 k events selected from the total sample 
of 667 k of registered and reconstructed events. This criterion essentially removes 
a contamination by interactions upstream of the target. For the event sample with 
the target removed, the selection reduces the number of events from 46 k to 17 k. 

4. Results 

This section presents results on inelastic and production cross sections as well as 
on differential spectra of 7r+ and tt~ mesons in p-l-C interactions at 31 GeV/c. 

4.1. Inelastic and production cross sections 

The total inelastic cross section is measured to b^^ 

cr„^; = 257.2 ± 1.9 ±8.9 mb . 

The production cross section was calculated from the inelastic cross section by 
subtracting the quasi-elastic contribution.'^l The result isP^ 

aprod = 229.3 ± 1.9 ± 9.0 mb . 

The production cross section is compared to previous measurements in Fig. |4j 

4.2. Spectra of 7r+ and tt~ mesons 

The TT"*" and 7r~ spectra presented in this section refer to pions produced in strong 
and electromagnetic processes in p+C interactions at 31 GeV/c. 
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Fig. 4. (Color online) Beam momentum dependence of the production cross section for p+C 
interactions. The NA61/SHINE result (filled square) is compared with previous measurements: 
Bellettini et al. (circle) p3 Carroll et al. (triangle JI3 and Denisov et al. (open squares) P^l For 
the NA61/SHINE point, the error bar indicates statistical and systematic uncertainties added in 
quadrature. The result from Ref.ESl was recalculated by subtracting from the measured inelastic 
cross section a quasi-elastic contribution at 20 GeV/c of (30.4 it 1.9 (sys)) mb. 



The spectra are presented as a function of particle momentum in ten intervals 
of the polar angle. Both quantities are calculated in the laboratory system. The 
chosen binning takes into account the available statistics of the 2007 data sample, 
detector acceptance and particle production kinematics. 

The final spectra are plotted in Figs. [5] and |6] For the purpose of a comparison 
of the data with model predictions the spectra were normalized to the mean tt^ 
multiplicity in all production interactions by dividing by <Tprod- This avoids uncer- 
tainties due to the different treatment of quasi-elastic interactions in models as well 
as problems due to the absence of predictions for inclusive cross sections. 



5. Conclusion 

This work presents inelastic and production cross sections as well as positively and 
negatively charged pion spectra in p-l-C interactions at 31 GeV/c on the 2007 thin 
(4% interaction length) carbon target. These data are essential for precise predic- 
tions of the neutrino flux for the T2K long baseline neutrino oscillation experiment 
in Japan. Furthermore, they provide important input to improve hadron production 
models needed for the interpretation of air showers initiated by ultra high energy 
cosmic particles. 

A much larger data set with both the thin and the T2K replica carbon targets 
was recorded in 2009 and 2010 and is presently being analysed. This will lead to 
results of higher precision for pions and extend the measurements to other hadron 
species such as charged kaons, protons, and A. Analysis of the data collected 
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Fig. 5. (Color online) Laboratory momentum distributions of 7r+ mesons produced in production 
p+C interactions at 31 GeV/c in different intervals of polar angle (8). The spectra are normalized 
to the mean Tr"*" multiplicity in all production p+C interactions. Error bars indicate statistical 
and systematic uncertainties added in quadrature. The overall uncertainty (2.3%) due to the 
normalization procedure is not shown. Predictions of hadron production models, FLUKA200di^ 
(solid line), UhQMdI. 3.5^ (dashed line) and VENUS4.li^ (dotted line) are also indicated. 
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p [GeV/c] p [GeV/c] 

Fig. 6. (Color online) Laboratory momentum distributions of tt~ mesons produced in production 
p+C interactions at 31 GeV/c in different intervals of polar angle (8). The spectra are normalized 
to the mean it~ multiplicity in all production p+C interactions. Error bars indicate statistical 
and systematic uncertainties added in quadrature. The overall uncertainty (2.3%) due to the 
normalization procedure is not shown. Predictions of hadron production models, FLUKA200diSl 
(solid line), UrQMd1.3.iE2I (dashed line) and VENUS4.ld2ni (dotted line) are also indicated. 
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with the T2K rephca target during the 2007 run is in progress!^ The new data 
will allow a further significant reduction of the uncertainties in the prediction of the 
neutrino flux in the T2K experiment. 
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